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ABSTRACT 

We have recently shown that both the Prugnicl-Simicn model and Scrsic's function (hereafter 
referred to as the Einasto model when applied to internal density profiles) describe simulated dark 
matter halos better than an NFW-like model with equal number of parameters. Here we provide 
analytical expressions for the logarithmic slopes of these models, and compare them with data 
from real galaxies. Depending on the Einasto parameters of the dark matter halo, one can expect 
an extrapolated, inner (0.01-1 kpc), logarithmic profile slope ranging from ~ —0.2 to ~ —1.5, 
with a typical value at 0.1 kpc around —0.7. Application of this (better fitting) model therefore 
alleviates some of the past disagreement with observations on this issue. We additionally provide 
useful expressions for the concentration and assorted scale radii: r s , r_2, r G , R c , r v - lr , and r max — 
the radius where the circular velocity profile has its maximum value. We also present the circular 
velocity profiles and the radial behavior of p(r)/a(r) 3 for both the Einasto and Prugniel-Simicn 
models, where a(r) is the velocity dispersion associated with the density profile p(r). We find 
this representation of the phase-space density profile to be well approximated by a power-law 
with slope slightly shallower than —2 near r = r_2- 

Subject headings: dark matter — galaxies: fundamental parameters — galaxies: halos galaxies: structure 
— methods: analytical 

1. Introduction 

1 Graham@mso.anu.edu.au 

In an interesting turn of events, Sersic's (1963, 
1968) 3-parameter function, developed to describe 
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the projected (on the plane of the sky) radial stel- 
lar distributions of galaxies, has been shown to 
also match the internal (3D) density profiles of 
simulated dark matter halos (Navarro et al. 2004). 
Intriguingly, this function was shown to provide a 
better fit than the 3-parameter NFW-like model 
with arbitrary inner power-law slope (Dicmand, 
Moore, & Stadel 2004; Merritt ct al. 2005, 2006). 
The functional form of Sersic's equation was inde- 
pendently developed by Einasto (1965, 1968, 1969) 
and used to describe the internal density profiles of 
galaxies (see also Einasto & Haud 1989 and Tenjes, 
Haud, & Einasto 1994). We shall therefore refer to 
this function as the Sersic model when applied to 
projected distributions, and as the Einasto model 
when applied to internal density profiles. 

One of the main concerns in (Merritt et al. 
2005, 2006) was whether the deprojected form 
of Sersic's function might provide a better de- 
scription of the density profiles than the Einasto 
model. Specifically, in Merritt ct al. (2006, here- 
after Paper I), the analytical approximation to the 
deprojection of Sersic's function given in Prug- 
niel & Simien (1997) was tested along with the 
Einasto model, the NFW-like model, and vari- 
ous other fitting functions. Overall, the Prugniel- 
Simien and Einasto models performed the best, 
providing a good description to both the galaxy- 
and cluster-sized halos built from hierarchical TV- 
body simulations, and also the monolithic cold col- 
lapse halos. 1 Curiously, while the Prugniel-Simien 
model provided the best fit to the cluster-sized ha- 
los, Einasto's model provided a better fit to the 
galaxy-sized halos. In this paper we explore some 
of the properties of these two models and some 
of the consequences they imply, and comparisons 
they enable, with real galaxies and galaxy clusters. 

In an effort to help clarify and unify the various 
parameters of the different models, Section 2 pro- 
vides relations between such quantities as effective 
radius R c , virial radius r v ; r , the radius where the 
logarithmic slope of the model equals —2, r_ 2 , and 
'concentration' as measured by observers and by 
modelers. 

In Section 3 we present the phase-space density 



1 In passing we note that the 3-paramctcr anisotropic 
Dchncn-McLaughlin (2005) model matched the galaxy- 
sized halos best, but it did not perform so well in describing 
the cluster-sized halos, and it was unable to describe the 
halos formed from spherical cold collapses. 



profiles, or more specifically, the density profiles 
divided by the cube of their associated velocity dis- 
persion profiles, showing how, for sufficiently large 
shape parameters n, both the Einasto model and 
the Prugniel-Simien model approximate a power- 
law <~ r~ 2 near r_ 2 . 

In Section 4 we derive the logarithmic slopes of 
the Einasto and Prugniel-Simien model, and com- 
pare these with real data. This is of particular 
interest because the innermost slope of these mod- 
els is considerably shallower than —1 and in fact 
equal to zero in the case of the Einasto model at 
r = 0. The inward extrapolation of these models, 
rather than the NFW-like model, therefore notice- 
ably reduces the disagreement between modelers 
and observers on this issue. 

Our findings are summarized in Section 5. 

2. The models: assorted radial scales and 
concentration 

We will discuss three (3-parameter) empirical 
models used for describing the internal density 
profiles of galaxies, clusters, and halos. Each 
model has three parameters and their application 
to simulated halos can be seen in Paper I. 

The first model is an adaptation of the Navarro, 
Frcnk, & White (1995, hereafter NFW) model to 
give a double power-law with an outer slope of —3 
and an arbitrary inner slope denoted by 7. The 
radial density profile, p(r), of this model can be 
written as 

P{r) = (r/r fl )T(l + r/r s )3-T (1) 

where r s is the scale radius at the density p Sl 
marking the center of the transition between the 
inner and outer power-laws with (extrapolated) 
slopes of —7 and —3. This function represents a 
restricted form of the more generic (a, f3, 7) model 
(Hcrnquist 1990, his equation 43; see also Zhao 
1996) and we shall refer to it as the (1, 3, 7) model. 
Setting (a,/3,7)=(l, 3, 1) yields the NFW model, 
while (1.5, 3, 1.5) gives the model in Moore et 
al. (1999). Because the (1,3,7) model has al- 
ready been studied in detail, our main focus will 
be on the models of Einasto (1965) and Prugnicl 
& Simien (1997). 

Einasto's model is given by the equation 

p(r) = Pc cxp{-d n [(r/re) 1 /" - l] } . (2) 
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The parameter n describes the shape of the den- 
sity profile. Larger values of n result in steeper 
inner profiles and shallower outer profiles. The 
quantity d n is defined to be a function of n such 
that p e is the density at the effective radius r 
which encloses a volume containing half of the 
total mass. A good approximation for n > 0.5 
is given by d n « 3n — 1/3 + 0.0079/n (Mamon 
2005, priv comm.), although we have used the 
exact value coming from L(3n) = 2 x j(3n, d n ), 
where 7 (xi, x<i) and T(x) are the incomplete and 
complete gamma functions, respectively (see Pa- 
per I). The mass profile (Mamon & Lokas 2005, 
their equation A2; Cardone et al. 2005, their equa- 
tion 11) is given by 



M(r) 



4irnrl p e e d " d n 



-3n 



7(3n, x). 



(3) 



We shall at times refer to the value of n from 
Einasto's model as nEin- 

The Prugnicl-Simien model can be expressed as 



p(r) = p' 



exp -b n (r/R e ) 



l/n 



with 



bn 



l(1 _ p) T(2n) 

2R e T(n(3 - p))' 



(4) 



(5) 



Once again, the parameter n describes the curva- 
ture of the density profile. The quantity b n is a 
function of n defined in such a way that R c is the 
effective radius containing half of the total mass 
when the 3D sphere defined by this density profile 
is seen in projection onto a 2D plane. Although 
we use the exact solution for b n , coming from 
r(2n) = 2x 7 (2n, b n ) (see Graham & Driver 2005), 
b n can be approximated by 2n — 1/3 + 0.009876/n 
for values of n > 0.5 (Prugniel & Simien 1997). In 
addition to n and R c , the third parameter which 
one solves for, p', is defined so that the volume- 
integrated mass from equation 4 is equal to the 
area-integrated mass of a Sersic function with the 
standard parameter set I C ,R C , and n (see Paper 
I)- 

The final term p is not a parameter but instead, 
like b n , another function of n. It is chosen to max- 
imize the agreement between the Prugnicl-Simien 
model and the deprojected Sersic model having 
the same parameters n and R c . A good match is 
obtained when p = 1.0 - 0.6097/n + 0.05463/n 2 



(Lima Neto et al. 1999; see also Paper I, their 
figure 13). The value of p is also responsible for 
determining the logarithmic profile slope at small 
radii. Setting p to zero, the Prugnicl-Simien model 
has the same functional form as Einasto's model. 

The internal density of the Prugnicl-Simien 
model at r = R c is given by p c = p'c~ bn , and 
the projected surface density at R — R c , de- 
noted by I e , can be solved for using equation 5. 
Comparisons of dark matter halos (fitted with the 
Prugnicl-Simien model) with real galaxies (fitted 
with Sersic's model) is obviously remarkably easy 
using this model. The associated mass profile 
(Lima Neto et al. 1999; Marquez et al. 2001) is 
given by the equation 

M(r) = 47rn J R e 3 pX" (3 " p)n 7 ((3 - p)n, Z) , (6) 

where Z = 6„(r/i? c ) 1 /". Expressions for the asso- 
ciated gravitational potential, force, and velocity 
dispersion can be found in Terzic & Graham 2005). 
In what follows, we shall at times refer to the value 
of n from the Prugnicl-Simien model as rips- 

In the following subsections we explore a num- 
ber of important radii associated with the above 
models, and address the issue of 'concentration'. 

2.1. The peak in the 4ttGV 2 p(r) [km s -1 ] 2 
profile, r 2 

As will be discussed in Section 4, the scale radii 
r c of the Einasto model, and the (projected) half- 
mass radii R c from the Prugnicl-Simien model, oc- 
cur where the logarithmic slope of the density pro- 
file is ~ —3. This can be quite far out, and so 
we define an additional radial scale. We do so by 
obtaining the radius where the profile 4irGr 2 p(r), 
which has units of velocity squared, has its max- 
imum. The integral of this profile gives the en- 
closed mass. 

For the (1, 3, 7) model (equation 1), this 
maximum occurs at a radius that we denote by 
r_ 2 .(i,3. 7 ), such that 



r-2,(i,3, 7 ) = ( 2 -7K 



7 < 2. (7) 



When 7 
dius r_ 



has r_ 
dius r_ 



2,(1,3,7) 
2,(1,3,7) 
2,(1,3,7) 



as in the NFW model 
r s , and when 7 = 



the ra- 

r s , and wnen 7 = 1.5 one 
— r s /2. It turns out that the ra- 
corresponds to the point where the 
logarithmic slope of the (1, 3, 7) density profile 
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equals —2, hence the adopted nomenclature. Sim- 
ilarly for the Einasto and Prugniel-Simicn density 
model, solving where the derivative of the profile 
4irGr 2 p(r) equals zero, one finds that these profiles 
also peak at the radius where their logarithmic 
slope equals —2. This is easy to understand when 
one notes that the solution to d[r 2 p(r)]/dr = 
leads to dp/p — —2dr/r. It is also easy to show 
that this corresponds to a maximum for any den- 
sity profile with a monotonically decreasing slope. 

For Einasto's model (equation 2), one has 



?*-2,Em 



2wEi; 

d n 



(8) 



and for the Prugnicl & Simien density profile 
(equation 4) one has 



r-2.PS 



nps(2 -p) 



(9) 



When nEin = 6, r_ 2 ,Ein ~ 0.10r o . When nps = 3, 
T-2, ps ~ 0.25i? c - (These representative values of 
the shape parameter have been taken from Pa- 
per I, which applied the above models to a number 
of simulated dark matter halos.) 

Evaluating Einasto's model (equation 2) at 
^-2,Ein to give the density p_2,Eim and express- 
ing r e in terms of r_ 2 ,Ein (equation 8), Einasto's 
model can be written as 2 

PEinM = (/9-2,Em)xexp{-2n[(r/r_2.Ein) 1/,l -l]}, 

(10) 

where 

P-2,Ein = PeC d - 2n . (11) 

This is the expression used in Navarro et al. 
(2004). 

Re-expressing Prugniel & Simien's model in 
terms of r_2, one has 3 



pps(r) = (P-2.ps) 

V-2.PS, 

x exp{-n(2 - p)[(r/r^ PS )^ n - 1]}, (12) 



where 



P-2,PS = P 



n(2 - p) 



exp{n(p-2)}. (13) 



2 For clarity, we have dropped the subscript 'Ein' from the 
parameter n. 

3 For clarity, we have dropped the subscript 'PS' from the 
parameter n. 



2.2. The peak in the circular velocity pro- 
file, 7* max 

While an isothermal density profile, p(r) oc r~ 2 , 
has a flat rotation curve, the radius where the log- 
arithmic slope of our (non-isothermal) density pro- 
files equals — 2, r_ 2 , does not coincide with the flat 
portion of the rotation curve, i.e. where the rota- 
tion curve has its maximum value. The rotation 
curves are simply given by the circular velocity 
profiles: u c irc 2 (= GM{r)/r), with M(r) defined 
previously. The maximum occurs at a radius r max 
which is larger than r_ 2 , and is shown in Figure 1 
as a function of both the effective radius and r_ 2 
for the Einasto and Prugniel-Simien models. 

The width near the peak of the circular velocity 
profile increases as the value of n increases. Den- 
sity profiles with larger values of n will approx- 
imate a flat rotation curve over a greater radial 
extent (in units of r_ 2 ) than profiles with smaller 
values of n (Figure lc and d). 

The radius r max can be obtained numerically 
by solving the expression for when the derivative 
of GM(r)/r equals zero. For the Einasto model, 
this amounts to solving 



7(3n,x) 



x 3n e- x /n, 



(14) 



with x = d n {r max /r e ) 1 / n . For the Prugnicl & 
Simien model, one needs to solve the expression 

7 (n(3 - p), Z) = Z n(3 -P^e- z /n, (15) 

with Z = &(r max /_R c ) 1 /™. The results are shown 
in Figure 2, with r max normalized against r_2- 

When n = 3.6 (and 2.9) in the Prugniel-Simicn 
model — the average profile shape for our galaxy- 
sized (and cluster-sized) CDM halos, see Paper I 
- f max ~ 2.17 r_ 2 (and 2.10 r_ 2 ). When n = 6 
(and 5.0) in the Einasto model, 2.21 r_ 2 

(and 2.16 r_2). This can be compared with the 
(7 = 1) NFW model for which r max is known to 
equal ~ 2.16 r_ 2 - 

2.3. Concentration and the virial radius, 

^*vir 

Given the rc-ncwed application of Einasto's 
model, which has the same functional form as 
Sersic's model that is used by observers in describ- 
ing projected distributions, it would seem relevant 
to inquire if we can also make use of the type of 
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Fig. 1. — Circular velocity profiles for a) the 
Einasto model (equation 2), and b) Prugniel & 
Simien's model (equation 4), for varying values 
of the profile shape n: n = 0.5 (solid lines), 
n = 1 (dotted), n = 2 (dashed), n — 4 (dash- 
dot), n = 10 (dash-triple dot). The lower panels 
show the same thing except the radius has been 
normalized against r_2 (the radius where the log- 
arithmic slope of the density profile equals —2) 
rather than the effective radii of the models. Panel 
c) Einasto's model, panel d) Prugniel & Simien's 
model, n — 0.5 (solid line); n = 1 (dotted line); 
n = 2 (dashed line); n — 4 (dot-dash line); n — 10 
(triple-dot-dash line). 
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Fig. 2. — The radius where the circular velocity 
profile peaks, divided by the radius where the log- 
arithmic slope of the density profile equals —2, is 
shown as a function of the density profile shape 
n for a) the Einasto profile (equation 2), and b) 
Prugniel & Simien's profile (equation 4). 



concentration indices that observers use. There 
are two flavors. 

The first is a ratio of two radii. While this may 
sound qualitatively similar to the NFW concen- 
tration, it is in fact fundamentally different. A 
classical example would be the ratio between the 
radii containing 50% and 25% of a galaxy's total 
light (Fraser 1972). In the case of a universal den- 
sity profile, with only a radial scale and a density 
scale, the ratio of radii containing 50% and 25% of 
the total (asymptotic) mass would always be ex- 
actly the same. That is, if we were to use the new 
concentration index r /r_ 2 , then if n is assumed to 
be constant (i.e., if a universal profile exists), this 
ratio will be the same for every profile (see equa- 
tions 8 and 9). A similar example comes from 
Butcher & Oemler (1984) who defined a galaxy 
cluster concentration index C — \og(R eo / R20) , 
where R x is the radius enclosing x% of the galaxies 
in a cluster. 

The second type of concentration index is a ra- 
tio of flux, compared to a ratio of radii, within two 
specified radii (Okamura, Kodaira, & Watanabc 
1984). An example is the flux within the radius 
containing half an object's total light divided by 
the flux within one third of this radius (Trujillo, 
Graham, & Caon 2001). But again, if the density 
profiles are universal, then such concentration in- 
dices will always have the same value. It is because 
real galaxies do not have universal light-profiles, 
i.e. a range of Sersic (i? 1 /™) indices are observed, 
that such concentrations indices work. 

Abraham et al. (1994) used a galaxy concentra- 
tion index defined as the ratio of flux within the 
radius containing an object's total light (rather 
than half its light) divided by the flux within 1/3 
of this radius. Now because, in general, galax- 
ies do not have well-defined edges but rather their 
light slowly peters out into the noise of the sky- 
background flux, deeper and deeper exposures 
yield increasingly larger total radii, and a flux ra- 
tio that tends to 1 for every galaxy. But because 
of the limited aperture sizes Abraham et al. (1994) 
used to define the total galaxy light, they obtained 
values different than 1. The quantity they mea- 
sured was thus a function of not only the light- 
profile shape, but how many R e they sampled in 
their largest aperture, as discussed in Graham, 
Trujillo, & Caon (2001). 

In a somewhat similar manner, the NFW con- 
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centration works because it too is dependent on 
the background noise, specifically, the mean mat- 
ter density of the universe. The virial radius, r v i r , 
is used to quantify the density of dark matter ha- 
los relative to the background. It is defined as the 
radius of a sphere containing an average matter 
density that is some specific number greater than 
the mean matter density in the universe. In Pa- 
per I we reported that our simulated halo profiles 
were computed using a value of 368. It is how- 
ever common to also see a value of ~337 when 
^baryon = 0.3, = 0.7, and h = 0.7 (e.g., Bryan 
& Norman 1998). Before the cosmological con- 
stant became fashionable, a value of ~178 was 
used for the flat Einstein de Sitter universe. 

Using equations 3 and 6 for the mass profile, one 
can (numerically) solve the following expression to 
obtain the virial radius in units of the scale radius 
r c and R c from the Einasto and Prugniel-Simien 
model, respectively. 

3M(r vir )/(47rr vir 3 ) = 337(p 

universe 

). (16) 

Halos do of course extend beyond r v i r (e.g., 
Maccio, Murante, & Bonometto 2003; Prada et al. 
2005). The results are shown in the top panels of 
Figure 3, where one can see how the virial radius 
is monotonically related to the density contrast 
between p e and (p U nivorso)- As r vlr increases be- 
yond the effective radius, the incomplete gamma 
function in the mass profiles starts to asymptote 
to a constant value and the slope in the figures 
tends to 1/3. 

From the relations connecting r c and R e with 
r_2 (equations 8 and 9), one can obtain the virial 
radius in units of r_ 2 . This is shown in the middle 
panels of Figure 3. The ratio r v i r /r_ 2 is referred 
to by modelers as the concentration parameter. It 
does not refer to the curvature or shape of the 
profiles, as observers might initially think, but is 
a measure of the density contrast of the halo rela- 
tive to the average background density of the uni- 
verse 4 . Obviously, if one did not wish to use the 
virial radius (see Maccio, Murante, & Bonometto 
2003 for an alternative prescription), then a sim- 
ilar 'concentration parameter' can be defined in 

4 Although these diagrams were created using an over- 
density factor of 337, this actual choice does not affect the 
(modeler's) concentration parameter's ability to act as a 
surrogate for the density scale p c or p—i- 



terms of /? e /(puniversc)- The slope at r = r v j r is 
shown in the lower panels of Figure 3 as a func- 
tion Of p c / (punivorso)- 

If one thinks of dark matter halos as icebergs, 
which can be lowered and raised relative to the 
background density of the universe, then profile 
universality means that one can use the offset be- 
tween (^universe) and either p e or p_ 2 or p' as a 
measure of 'concentration'. But if a range of pro- 
file shapes exists, i.e. different n (a in the notation 
of Navarro et al. 2004), the difference between p e 
and p-2 and p' will depend on the profile shape. 
What this means is that the concentration one 
measures will depend on where one samples the 
halo's density. This is important because the halo 
density is thought to reflect the mean density of 
the universe when the halo formed, and is thus a 
measure of the collapse redshift of the halo. 

It would be of value to explore whether or not 
the use of p e and M(r e ) (and M(R C ) in the case 
of the Prugniel-Simien model), rather than p_ 2 
and M v i r may account for some of the scatter in 
diagrams plotting concentration versus halo mass, 
or equivalently, scale radius versus scale density. 

Finally, we note that the use of a Petrosian- 
style radius (Petrosian 1976; Graham et al. 2005), 
such that the mean density inside of some radius 
divided by the density at that radius equals some 
constant value, is not suitable in the case of struc- 
tural homology. This is because such a radius 
will be equal to the same fractional number of 
scale radii (r c or r_ 2 ) for every halo. That is, 
a Petrosian-like radius will just be a re-expression 
of the scale radius. 

3. On the power-law nature of p/a 3 

There has been recent interest in the the 
pseudo phase-space density profiles represented 
by p(r)/cr(r) 3 . This quantity appears to be well 
approximated by a power-law r~ a , with a w 1.94 
(Taylor & Navarro 2001; Ascasibar et al. 2004; Ra- 
sia, Tormen, & Moscardini 2004; Sota et al. 2006; 
Barnes et al. 2006). Independently of any model, 
Dchncn & McLaughlin (2005) found a best-fit 
value of a = 1.92 ± 0.01 using the halos A09-F09 
and G00-G03, which we studied in Paper I. 

Figure 4 shows the ratio p(r)/o-(r) 3 for the 
Prugniel-Simien density profile (equation 4) cou- 
pled with its spatial (i.e., not projected) velocity 
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Fig. 3. — Panels a) and b) show the virial radius (normalized against the effective radius r and R e , respec- 
tively), as a function of the scale-density p c relative to the average background density of the universe for 
Einasto's model (equation 2) and the Prugniel-Simien model (equation 4), respectively. Panels c) and d) are 
similar to panels a) and b) but with the virial radius normalized by the radius where the logarithmic slope 
of the density profile equals —2. Panels e) and f) show the associated, negative, logarithmic slope of the 
density profile at the virial radius, denoted by 7 v i r - n — 0.5 (solid line); n = 1 (dotted line); n = 2 (dashed 
line); n = 4 (dot-dash line); n = 10 (triple-dot-dash line). 
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dispersion profile given in Terzic & Graham (2005, 
their equation 28). Isotropy in velocity space has 
been assumed. As can be seen, the profiles are not 
exactly featureless power-laws, but for n > 4 the 
departure from a power-law, over the radial range 
shown, is less than about 20%. 

For the Einasto density profile (equation 10), 
the spatial velocity dispersion profile can be ob- 
tained by integrating the isotropic Jeans equation 
of hydrostatic equilibrium 



<y s 2 {r) 



G 
p(r) 



(17) 



Expressing r c in terms of r_ 2 (equation 8) in the 
mass profile M(r) (equation 3), one obtains 



<T s 2 (x) 



GM tot (2n) 



1+71 



r-2 



2 r(3n) 

POO 

/ S-"- 1 e- s 7 (3n,x)dx, (18) 

J X 



where x — d n (f /r c ) 1/>n — 2n(f/r_ 2 ) 1 ' / ". Integra- 
tion to infinity for this expression 5 is avoided by 
making the substitution x = {xj cos 9), such that 
dx/d8 = xsinO/ cos 2 0, giving 



er s 2 (x) 
t/2 



GM t , 



(2n) 1+n e x x~ 



r-2 



r(3n) 



x / tantf cos™ ^ e" x/cose 7 (sn, d9.(l9) 
J Q \ cos 9/ 

Figure 5 shows p/cr s 3 for the Einasto density pro- 
file. 

From the residual profiles in Figures 4 and 5, 
one can see that, over the radial range 0.1 < 
r/r_2 < 10, a slightly shallower slope than —2 
exists for 4 < n < 10 (Einasto case, Figure 5) and 
2 < n < 10 (Prugniel-Simien case, Figure 4) 6 . Al- 
though at r = r_ 2 , — dlogp/dlogr = 7 = 2, for 
«Ein = 6,7 Ein (0.1r_ 2 ) = 1-36 and 7Ein(10r_ 2 ) = 
2.94, and for n PS = 3, 7 PS (0.1r_ 2 ) = 1-62 and 



5 Equation 23 in Cardone et al. (2005) for the velocity dis- 
persion profile contains a typo such that the term in the 
exponential should be (+2/j)(r /r-2) 1 rather than (—2/7), 
where their 7 equals 1/n. Their figure 5 is however correct. 

6 Had the slope been exactly —2 in Figures 4 and 5, it would 
have implied an outer logarithmic slope in the density pro- 
file that decayed in an oscillary manner about a value of 
—2 (Dchnen & McLaughlin 2005, their figure 2), at odds 
with the data in Figure 3 from Paper I 
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Fig. 4. — Prugniel-Simien density profile p(r) 
(equation 12) divided by the cube of its spa- 
tial (i.e., not projected) velocity dispersion pro- 
file cr s (r) (Terzic & Graham 2005, their equa- 
tion 28). The curves are such that: n = 0.5 
(solid line), n = 1 (dotted), n — 2 (short-dashed), 
n = 4 (dash-dot), n = 10 (long-dash). The curves 
asymptotically approach a line having slope —2 
(shown by the solid straight line) as n — > 00. One 
obtains the same asymptotic behavior using the 
projected velocity dispersion profile. The lower 
panel shows the difference between the curved pro- 
files and the line of slope —2, divided by the den- 
sity of the curved profiles. 
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7ps(10r_2) = 2.56. The slope of the density pro- 
file does therefore change over this radial range. 
This can be appreciated in Figures 6 and 7 which 
show the negative, logarithmic slope of the (cube 
of the) velocity dispersion profile, the density pro- 
file, and p/o- s 3 • For large values of n, p/o- s 3 ~ r~ 2 . 

4. Density profile slopes 

De Blok (2005) has argued that the inner pro- 
file slopes of simulated, dark matter halos are in- 
consistent with observations of dark matter domi- 
nated galaxies (Moore et al. 1999; Salucci & Burk- 
crt 2000; Marchesini ct al. 2002; de Blok, Bosma, 
& McGaugh 2003; Gentile et al. 2006; Goerdt et 
al. 2006). He reports that the inner density pro- 
files of low surface brightness (LSB) galaxies have 
logarithmic slopes significantly shallower than — f 
at a radius of 0.4 kpc. This is important because 
it suggests a possible problem with hierarchical 
ACDM simulations of dark matter halos, which, 
at least from (1, 3, 7) model fits, typically have 
inner slopes steeper than —1. 

While there is presently no consensus as to 
why such a disagreement exists, some of the ap- 
parent discrepancy may arise from either bary- 
onic processes which modify the dark matter pro- 
file (e.g., Mashchenko, Couchman, & Wadslcy 
2006), or from systematic biases in measuring in- 
ner slopes from HI and Ha long-slit observations 
(van den Bosch et al. 2000; Swaters et al. 2003a; 
Spekkens, Giovanelli, & Haynes 2005, but see de 
Blok 2003). For example, non-circular motions 
can make galaxies appear less cuspy than they re- 
ally are. Significant, in the sense of non-zero, non- 
circular motions are indeed present in many LSB 
galaxies where high resolution 2D velocity fields 
are available (e.g., Swaters et al. 2003b; Blais- 
Ouellette et al. 2004; Coccato et al. 2004; Simon 
et al. 2005). However, on their own, these do not 
explain the observed difference in slope (de Blok 
ct al. 2003; Gentile et al. 2006). But in combi- 
nation with gas pressure support and projection 
effects, Valenzuela et al. (2005) argue that this 
may account for the relatively shallow slopes in 
observations. Using extensive simulations of ob- 
serving and data processing techniques, Spekkens 
et al. (2005) also report how measurements from 
long-slit, optical spectra of halos with inner slopes 
of —1 can result in "observed" slopes consistent 



with values around —0.25 ± 0.15. The apparent 
success of the flat-core Burkert (1995) model may 
then be an artifact of observational biases. Higher- 
resolution gamma-ray studies of dark-matter dom- 
inated galaxies may, in the future, help to resolve 
the current cusp-core controversy (e.g., Profumo 
& Kamionkowski 2006; Lavalle et al. 2006). 

The inward extrapolation of simulated density 
profiles using empirical models which have steep 
(asymptotic) inner power-laws may also be partly 
responsible for the mismatch (see., e.g., Kravtsov 
ct al. 1998). As noted by Navarro et al. (2004) 
and Stoehr (2006), empirical models with shallow 
inner slopes, such as Einasto's model, not only 
match the simulated data down to 0.01 r v j r but 
could potentially resolve the apparent dilemma at 
smaller radii (but see Diemand et al. 2005, who 
find a slope of —1.25 at 0.001 r v ; r in a highly- 
resolved, cluster-sized halo). That is, ACDM cos- 
mology and the various A^-body simulations them- 
selves may in fact be fine, but the empirical models 
used to parameterize the CDM halos may fail at 
small radii. 

Here we examine the slope of the various empir- 
ical models, and compare these with observations 
of real, dark matter dominated galaxies. 

In the case of the (1, 3, 7) model (equation I), 
the slope is given by 

7(i,., 7 ,(r)-^^=7+(3-7)/(l + r./r). 



dlogr 

For small values of r/r s , 
7(1,3,7) ( r ) ~7+ (3 -7); 



(20) 



r«r s , (21) 



which, as expected, asymptotically approaches 7 
as r — > 0. Figure 8 shows the negative, logarithmic 
slope as a function of radius for a sample of (1, 3, 
7) models with 7 =0.5, 1.0, and 1.5. One can see 
that the negative logarithmic slope of the profiles 
are practically equal to 7 at r < 0.01r s . What 
should also be realized is that, although the (1, 
3, 7) models do have continuously curving slopes 
from 0.01 to 1 r v ; r , they don't have the correct con- 
tinuously curving slope to match the CDM halos 
as well as Einasto's model can or the Prugnicl- 
Simien model can. 

The negative logarithmic slope of Einasto's 
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model (equation 2) is given by 



7Ein(r) 



loge 



d[logp(r)] d n ( r 



dr 



n 



l/n 



(22) 

which is approximately 3(r/r ) 1 /™ for n > 1 (see 
Figure 11 in Paper I). One can also see that when 
r = r c , the negative logarithmic slope of the den- 
sity profile is approximately 3. From Figure 3 of 
Paper I it is clear that r c will occur at a large ra- 
dius. For fixed values of n, Figure 9a shows how 
the negative logarithmic slope decreases monoton- 
ically as the radius r decreases. 

The negative logarithmic slope of Prugniel & 
Simien's model (equation 4) is given by 



7Ps(r) 



K 
n 



l/r, 



+p- 



(23) 



From b n w 2n - 1/3 + 0.009876/n for n > 0.5 
(Prugniel & Simien 1997), 7P s(r) w 2{r/R c ) 1 / n + 
p, and thus ~ (2 + p) at the effective radius R e . 
For large n, p — > 1 and 7ps(?"e) - * 3, as is the 
case with 7Ein(''e)- These profile slopes are shown 
in Figure 9b as a function of radius for different 
values of the profile shape n. 

Using equations 8 and 9, one can reformulate 
the above equations to obtain 



7Ein(r) = 2(r/r_ 2 ,Ei„) 1/ ", 



and 



7P S (r) = (2-p) 



r~2,ps 



l/n 



+ P- 



(24) 



(25) 



As r — > 0, 7Ein — » 0, apparently in fair agree- 
ment with the observations of real galaxies re- 
ported in, for example, Simon et al. (2003) and 
de Blok et al. (2003), who find a negative loga- 
rithmic slope of 0.2 ± 0.2 (but see Section 4.1). In 
the case of the Prugnicl-Simicn model, as r — > 0, 
7ps -» P (= 1-0 - 0.6097/n + 0.05463/n 2 ). Re- 
sults from Paper I gave galaxy (and cluster) pro- 
file shapes ranging from nps = 3.14 to 4.55 (and 
from nps = 2.19 to 3.47), suggesting a range of 
central (r = 0), negative logarithmic profile slopes 
for the Prugniel-Simicn model of 0.81-0.87 (and 
0.73-0.83). These slopes are considerably shal- 
lower than the mean (± standard deviation) value 
7 = 1.32 ± 0.19 (and 1.15 ± 0.16) obtained from 



the (1, 3, 7) model fits (Paper I). They are also 
in excellent agreement with theoretical expecta- 
tions based on phase-space arguments which sug- 
gest that CDM density profiles should have central 
cusp slopes equal to 0.75 (Taylor & Navarro 2001, 
see also Hansen & Stadel 2005 and An & Evans 
2006). 
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Fig. 5. — Same as figure 4 except that Einasto's 
density profile (equation 10) has been used here. 



icr 2 io - ' 1 10 

r/r_ 2 



Fig. 6. — Negative, logarithmic slopes associated 
with the Prugniel-Simien model (equation 12). 
The curves in the lower panel equal the curves in 
the middle panel minus those in the upper panel. 
See section 3 for details. 
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Fig. 7. — Same as figure 6 except that Einasto's 
density profile (equation 10) has been used here. 



o 

T 



log(r/r g ) 

Fig. 8. — Negative logarithmic slope of the (1,3, 
7) model (equation 1) when 7=0.5, 1.0, and 1.5 
(bottom, middle, and top curves, respectively). 
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log(r/R) 
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Fig. 9. — Panel a) Negative logarithmic slope of 
the Einasto r 1 /" profile (equation 22) as a func- 
tion of normalized radius r/r c for different values 
of n = 0.5, 1, 2, 4, and 10. Panel b) Negative loga- 
rithmic slope of Prugniel & Simien's density profile 
(equation 23) MS cX function of normalized radius 
r/R e . Panels c) and d) are the same as a) and b) 
except that the radius has now been normalized at 
r -2,Ein an d f*-2.PSi respectively (see equations 8 
and 9). As r — > 0, 7Ei n — > while 7ps — > p. 
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4.1. Slope comparison with real galaxies 

For a more meaningful comparison between 
model halos and observations of real galaxies, ob- 
servers and modelers should report profile slopes 
as a function of radius and perform their compar- 
isons at the same radii. Remarks in the literature 
that higher resolution rotation curves tend to show 
the greatest departure from an inner logarithmic 
slope of —1 (or —1.5) are somewhat beguiling. Be- 
cause such measurements of the inner profile slope 
in real galaxies were often made at radii smaller 
than those typically probed by ACDM simula- 
tions, they do not provide a particularly strong 
constraint or check on the simulations. They do 
however provide a check on any empirical fitting 
functions whose inward extrapolation does not fall 
below some fixed slope, such as — 1. Addressing 
this issue, de Blok et al. (2005) has compared real 
and simulated systems at 0.4 kpc. He found that 
the density profiles implied by the best-fitting, 3- 
parameter function used by Hayashi et al. (2004; 
equation 8 from Rix et al. 1997) to model the ve- 
locity profiles of LSB galaxies have slopes which 
are inconsistent with a value as steep as —1, and 
thus also with the average value of ~ —1.2 that is 
typically reported for simulated halos. 

For a mean value of riEin ~ 6 (from the Einasto 
r 1 /" fits in Paper I), a negative logarithmic slope 
of 0.5 occurs at 2.4 x 10~ 5 r o (2.4 x 10~ 4 r_ 2 )- This 
is about 10 pc for a galaxy halo with R c = 400 kpc, 
and corresponds to 0.12 arcsec at the distance of 
the Virgo cluster (17 Mpc). At 0.1 kpc, a typi- 
cal value at which observers measure the slope of 
the mass-density profile in real galaxies (see Fig- 
ure 10), one would expect to find a negative loga- 
rithmic slope equal to 0.73 for this halo; in perfect 
agreement with the mean slope obtained by Simon 
et al. (2005) for a sample of real galaxies. At 0.4 
kpc, one has 7Ei n ~ 0.92. If n = 5 and R c = 200 
kpc, then at 0.4 kpc one has 7Ei n = 0.85, and at 
0.1 kpc 7Ein = 0.64, consistent with the data from 
Swaters et al. (2003a). 

Figure 10 shows the innermost, resolved, loga- 
rithmic slope from the density profiles (assuming 
a minimum stellar disk) of 70 faint, LSB galaxies 
thought to be dark matter dominated (de Blok & 
McGaugh 1997; but see the warning 7 in Graham 

Although many are, not all LSB galaxies are particularly 
dark-matter dominated. For example, UGC 3137 and 



2002), plotted against the physical radius at which 
the slope was measured, i?mncr- This figure has 
been adapted from de Blok (2003, his figure 3) 8 . 
In order to compare how well the new density mod- 
els perform, it is necessary to plot several profiles 
with differing scale radii — which amounts to a 
horizontal shift of the curves in Figure 10. 

Einasto's r 1 /" model appears capable of match- 
ing the data reasonably well, depending on the 
combination of scale radius and profile shape n. 
However, for the halos studied in Paper I, bounded 
by the curves shown in Figure 10, the best-fitting 
Einasto models do not have negative, logarith- 
mic slopes shallower than ~0.4 at radii >0.1 kpc. 
This is at odds with roughly half of the galax- 
ies from de Blok et al. (2001) and de Blok & 
Bosma (2002), but largely in agreement with the 
data from Swaters et al. (2003a) and Verheijen 
(1997). Clearly, the apparent inconsistency be- 
tween the inner profile slope of dark matter halos 
generated from ACDM A^-body simulations and 
observations of real galaxies is reduced upon re- 
placement of the NFW model with the (better fit- 
ting) Einasto model. What is also apparent is that 
one can expect a range of different slopes inside of 
1 kpc; although this in itself does not imply a non- 
universal density profile. The largest study to date 
of (165) low-mass galaxies found inner logarithmic 
slopes ranging from -0.22 ± 0.08 to -0.28 ± 0.06 
for various subsamples of the data (Spekkens et al. 
2005). However, after extensive testing, these au- 
thors concluded that, due to biases in the analysis 
of long-slit spectra, the data is in fact consistent 
with inner logarithmic slopes ranging from to 
-1. 

The extrapolation of the Prugniel-Simicn model 



UGC 5750 from de Blok & Bosma (2002) have faint, central 
B-band surface brightnesses of 24.1 and 23.5 mag arcsec -2 
respectively, yet their total mass (stars, gas, dark matter) 
within 4 scale-lengths (=v IQ t 2 ih / G) divided by their flux 
within this radius (equal to 91% of the total, exponential, 
disk flux) gives a solar Mtot/Ln ratio of only 13 and 11, 
respectively. Typical M^/Lg ratios for Sa-Sd galaxies, 
within ~4 scale-lengths, are 3 to 7 (Roberts & Haynes 
1994). Baryonic processes (Weinberg & Katz 2002) might 
therefore be important here, especially if fractionally more 
HI gas exists in LSB galaxies. 

Figure 10 differs slightly from figure 3 in de Blok (2003) 
because we have included all 15 data points from Swaters 
et al. (2003a), and we have correctly reversed the symbols 
used to differentiate the data from de Blok et al. (2001) 
with that from de Blok & Bosma (2002). 
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inside of ~1 kpc — the extent to which our simu- 
lations (Paper I) provide meaningful data — does 
not do so well at matching the observations of real 
galaxies (Figure 10b). This remark does however 
overlook the previously noted analysis of Spckkcns 
ct al. (2005). In passing we note that Demarco et 
al.'s (2003) analaysis of 24 galaxy clusters using 
the Prugnicl-Simien model yielded a mean inner 
slope of -0.92. 

Although not explored in this work, the power- 
law exponent p in the Prugniel-Simien model can 
be modified. As we saw in Figure 9, as r — > 0, the 
negative logarithmic slope of this model tends to 
the value p. If one was to reduce the value of p, 
then one would acquire shallower inner slopes. As 
we noted before, if one reduced the value of p to 
zero, then one would obtain Einasto's model. 

Ideally, rather than simply plotting the inner 
profile slope versus the radius in kpc at which 
the slope has been measured, one should factor 
in that galaxies possess a range of sizes, i.e. scale 
radii. For example, 0.4 kpc may correspond to 0.5 
scale radii or 0.05 scale radii. Although neither 
the Einasto nor Prugniel-Simien models have yet 
been applied to the rotation curve data of the 70 
galaxies shown in Figure 10, a pseudo-isothermal 
model has been fit to most of these galaxies. This 
simple model can be written as 

*> " TTW' <26) 

where r c is the core radius and po the central den- 
sity. Figure 11 shows the logarithmic slope of 
this model, together with data from de Blok et 
al. (2001), de Blok & Bosma (2002), and Swaters 
et al. (2003a) for which values of the scale radius 
r c were available. The scattering of points, rather 
than following the curve, suggest that the data 
do not behave according to the pseudo-isothermal 
model, and/or an underestimation of i?i nne r rel- 
evant to where the slope was measured. Note 
though that the pseudo-isothermal model is an ex- 
treme model with a somewhat large, flat inner den- 
sity profile: models based on recent observations 
favor a slightly steeper slope of — 0.2±0.2 (e.g. de 
Blok ct al. 2001, 2003), while others find a slope 
scattered around —0.73 ± 0.44 (e.g. Simon et al. 
2005; Swaters et al. 2003a). Steeper cusp mod- 
els would provide a better fit to the data shown 
in Figure 11. It would be of interest to obtain the 





a) Einasto " 




J b) PS97 " 






'Wim t 1 







l°g(JWAP C ) lo g( R innerAP C ) 

Fig. 10. — Adaptation (see footnote 8) of figure 3 
from de Blok (2003). Data points show the log- 
arithmic slope of the density profile (assuming a 
minimum stellar disk) at the innermost resolved 
radius for a sample of 70 real galaxies. Open cir- 
cles: Verheijen (1997); filled circles: de Blok ct al. 
(2001); filled squares: de Blok & Bosma (2002); 
stars: Swaters ct al. (2003a). Over-plotted in 
the left panel are the profile slopes from Einasto's 
r 1 /" model (equation 10) with r_ 2 = 10 kpc (dot- 
ted lines) and 100 kpc (solid lines) for n=4 (up- 
per curve), 6 (middle curve), and 8 (lower curve). 
The right panel shows the same thing but for the 
Prugniel-Simien model (equation 12) with n=2 
(upper curve), 3 (middle curve), and 4 (lower 
curve). 
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best-fitting Einasto radii r e and profile shapes nEin 
for these galaxies, which would allow one to com- 
pare how well the observed inner slopes correlate 

With -RinnerAc- 



5. Summary 

We have provided expressions to relate the half- 
mass radii of the Einasto and Prugniel-Simicn 
models to a) the radius, r_2, where the logarith- 
mic slope of the density profile equals — 2, b) the 
virial radius, r v ; r , and c) the radius where the as- 
sociated circular velocity profile has its maximum 
value, r max . 

We have shown the dependence of the 'concen- 
tration' terms r v i r /r_2 and r v - lx /r c (and r vlr /R e ) 
on the ratio p e /(p U niv), where p e and r e (and R e ) 
are the Einasto (and Prugniel-Simien) scale den- 
sity and half-mass radius. We also show how the 
slope of these models at r v ii depends solely on the 
ratio p c /(Puniv) (Figure 3). 

Over the radial range 10~ 2 < r/r_2 < 
4 x 10 1 , we have shown both the Einasto and 
Prugniel-Simien models possess the property that 
p(r)/a(r) 3 can be roughly described by a power- 
law r~ a with the value of a slightly less than 2 for 
profile shapes n equal to or greater than the best- 
fitting values reported both here and elsewhere. 

Analytical expressions for the logarithmic slope 
of the Einasto and Prugniel-Simien models have 
been derived, and the slope expected from the 
inward extrapolation of these models, inside of 
~0.01 r v j r , is compared with that from observa- 
tions of real galaxies. The innermost (r = 0) 
slope of the Prugniel-Simicn model (0.73-0.87), as 
currently defined, appears too steep to match all 
the galaxy data, but agrees with theoretical ex- 
pectations for a slope of —0.75 (Taylor & Navarro 
2001) and -0.78 (Austin et al. 2005). Future work 
should explore the optimal value of the quantity p, 
the inner logarithmic profile slope in the Prugniel- 
Simicn model. Setting p = 0, one recovers the 
Einasto model, which appears capable of matching 
the inner profile slopes observed in real galaxies 
(Figure 10). Indeed, the typical value of ~ —0.7 
at 0.1 kpc in our CDM halos agrees well with the 
galaxy data from Swaters et al. (2003a) and Si- 
mon et al. (2005), but is steeper than the value 
-0.2 ± 0.2 reported by de Blok et al. (2003) and 
others. We also note that, at present, the pseudo- 
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Fig. 11. — Solid line: logarithmic slope of the 
pseudo-isothermal model (equation 26). The data 
points are from observations of real galaxies (as- 
suming a minimum stellar disk). Open circles: 
Verheijen (1997); filled circles: de Blok et al. 
(2001); filled squares: de Blok & Bosma (2002); 
stars: Swaters et al. (2003a). If these galaxies 
were described by the pseudo-isothermal model, 
they would follow the curve. 
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isothermal model appears inconsistent with the 
galaxy data (Figure 11). 
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